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Abstract 

Based on simple physical and geometric assumptions, we have calculated the mean surface molecular 
density of spiral galaxies at the threshold between star formation induced by cloud-cloud collision and 
spontaneous gravitational collapse. The calculated threshold is approximately logS cr i t ~ 2.5, where £ M Q • 
pc~ 2 is the observed surface mass density of an assumed flat gas disk. Above this limit, the rate of molecular 
cloud collisions dominates over spontaneous molecular cloud collapse. This model may explain the apparent 
discontinuity in the Schmidt law found recently at 2 < logS < 3. 
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1. Spontaneous and collisional star formation 

The onset of star formation within an ensemble of 
molecular clouds is thought to have two modes, spon- 
taneous (stochastic) gravitational collapse (Bonncll ct 
al. 2004; Krumholz & McKee 2005), and collisions be- 
tween molecular clouds (e.g., Tan 2000). Most studies to 
date have focused on one of these two possibilities. In 
actual physical situations which cover a range of physical 
parameters, these two processes should both have effect 
on the overall star formation in galaxy disks. In regions 
of higher molecular density, collisions should become dom- 
inant over spontaneous formation, because collisions will 
occur on timescales shorter than spontaneous star forma- 
tion. In regions with scarse cloud number density, col- 
lisions should also be rare and therefore star formation 
from other stochastic means should dominate. The tran- 
sition from spontaneous to collisional star formation may 
have observed effects on the relation between molecular 
gas content and star formation rate (SFR) . This expected 
transition has not been clearly observed yet, however. 

There is growing observational evidence that cloud col- 
lision can trigger star formation (Loren 1976; Scoville 
et al. 1986; Hasegawa et al. 1994; Koda & Sofue 2006). 
These studies show that cloud collisions are seen mainly 
at high density regions, and therefore merging and inter- 
acting galaxies are candidates for galaxies forming stars 
mainly by collision. A comparison of these galaxies 
with normal galaxies have been conducted by Young et 
al. (1986). They have found that merging/interacting 
galaxies exhibit SFRs that are an order of magnitude 
higher than normal galaxies of the same molecular gas 
density. Similarly, central regions and spiral arms of nor- 
mal galaxies, which generally have high molecular den- 
sities, may include star formation induced by collision. 
However, quantitative calculations which give direct (ob- 
servable) parameters on exactly where in parameter space 



star formation transits from spontaneous to collisional 
mechanisms have not been proposed. 

2. Model 

We consider a spiral galaxy that can be modelled as a 
thin disk. We take a cylinder of height 2d and base area 
S within the same plane as the galaxy disk (see schematic 
figure 1). We can think of d as the scale height, and S 
as the area of the deprojected beam of the observation. 
Within this cyclinder are n spherical molecular clouds 
with diameter D, mass M and mean molecular density 
p. Assuming that the observer can see through the disk 
and derive a surface molecular gas density £, the obvious 
equations are: 

nM , . 

S = ^ (1) 
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These molecular clouds are moving about the cylinder 
with an inter-cloud velocity dispersion a. 

For further analysis, we have several simple assump- 
tions: 

• star formation occurs on the free fall timescale, de- 
pendent only on the initial density 

• molecular clouds are moving about the thin cylinder 
in random motion 

• star formation in galactic disks is a two dimensional 
process. 

The third assumption is motivated from the fact that 
giant molecular clouds have size comparable to the scale 
height d (Stark & Lee 2005; Malhotra 1994). Wc can 
therefore assume that all processes occur in two dimen- 
sions and we can hereafter ignore d hereafter. 
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Table 1. Derived values of £ c rit 
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Fig. 1. Schematic picture of model. The base of the cylinder 
is parallel to the galactic disk. Each of the molecular clouds 
have size D, mass M, density p, and velocity a. 

We have 

/ 3tt \* 
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where iff is the free fall timcscale and G the gravitational 
constant. Collisions between molecular clouds should have 
a mean free path A m f p , for two dimensions, given by 

Amf ^7lb (4) 

where N = n/S is the number density of molecular clouds. 
A GMC harboring spontaneous star formation will move 
a linear distance of atg during its formation of stars. In 
a regime where spontaneous star formation is dominant 
over collision, parameters should satisfy the relation 

atg <C A m fp 

so that the GMC does not collide with other GMCs during 
formation. On the other hand, a regime where collision is 
dominant should satisfy 

atg ^ A m f p . 

Therefore, if a threshold exists between these two regimes, 
parameters should satisfy the equation 

crt ff ?s A m f p . (5) 

Substituting equations (1) through (4) into (5) and denot- 
ing the critical £ as S C rit, we obtain 
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Now we apply actual values to equation (6). For typ- 
ical Galactic molecular clouds with D ~ 60 [pc] and p ~ 
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Notes- D, the diameter of molecular clouds, are taken 
so that it may represent smallest (lOpc) and largest 
(300pc) clouds, a from literature ranges from 4 km/s 
(Liszt k Burton 1981; Clemens 1985) to 7km/s (Stark & 
Brand 1989), but we take values from 1 to 10 km/s to 
demonstrate the range of values £ cr it may take. 
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and inter-cloud velocity disper- 
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sion er^4[km-s~ 1 ], we obtain log£ cri t = 1.91 [M Q pc~ 
For observed molecular surface densities above this value, 
collisions dominate. 

We can further reduce the number of variables. Dame 
et al. (1986) have found that for Galactic molecular 
clouds, p scales with D such that logp = 2.80— 1.321og£>. 
Substituting this relation into equation (6) gives 
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(7) 



It is apparent from equation (7), that the critical density 
is practically independent of GMC diameter. The critical 
surface density is dependent mainly on cloud velocity dis- 
persion a, which is generally taken to be below 10 km-s . 
Table 1 lists values of £ C rit for a range of parameters to 
evaluate the lowest and highest cases. 

Errors for this value can be evaluated based on the un- 
certainty of the parameters, a derived from other stud- 
ies, numerical and observational, have agreed-upon val- 
ues between 1 [km/s] and 10 [km/s]. The value of D, 
which must be taken so that D is representative of the 
size which dominates star formation, is between 10 [pc] 
and severalx 100[pc], although uncertainties in D have lit- 
tle effect on E cr it- On the average, we may safely apply 
errors so that S C rit = 2. 
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4. Comparison with observation 

The empirical relation between gas density and SFR, 
often referred to as the Schmidt law, has long been stud- 
ied, and the relation of the form Ssfr oc E w , where Esfr 
is the area averaged SFR, has been extensively studied. 
Details can be found elsewhere (e.g., Buat et al. 1989, 
Buat 1992, Kennicutt 1998, Rownd and Young 1999,Wong 
and Blitz 2002, Boissier et al. 2003, Yao et al. 2003, Gao 
and Solomon 2004, Heyer et al. 2004). This law, however, 
owes its dynamic range to the contrast between starburst 
galaxies and normal spiral galaxies. The starbursts have 
higher molecular density whereas normal spirals have con- 
siderably lower molecular density, so the combination of 
these two types of galaxies contributes dominantly to the 



Usage of pas j 00. els 



3 



>, 

Q_ 

-52 
z; 

C/3 



V) 

W 

□) 
o 




-1 01234567 
log E (Msun/sqpc) 



Fig. 2. The Schmidt law, 
al. (2005). Abscissa is log £ 
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galaxies, and open squares represent circumnuclear star- 
bursts. The center vertical line is the value for £ cr it = 2.8, 
where the dashed lines represent errors as stated in text. 

linear fit of the Schmidt type power law. 

Any systematic differences in the Schmidt law between 
these two types of galaxies are ignored by fitting both with 
the same function. Indeed, a possible significant difference 
in these two types of galaxies were found by Komugi et 
al. (2005), by the study of high molecular density central 
regions of normal spirals. The Schmidt law, with nor- 
mal galaxies and circumnuclear starbursts superposed, is 
shown in figure 2. 

An apparent discontinuity exists in the two sequences, 
and the critical surface density E C rit = 2. 8+j' J , shown as a 
vertical line, coincides with the transition. Assuming that 
this apparent transition is real, it may be explained by 
change in the physical nature of star formation, namely 
the transition from spontaneous to cloud-cloud collisional 
star formation, induced by the increase in molecular den- 
sity. 

It must be pointed out that the SFR for the circumnu- 
clear starbursts are derived from far-infrared (FIR) lumi- 
nosities, and those of normal galaxies from internal extinc- 
tion corrected Ha luminosities (see Komugi et al. (2005) 
for details). 

These two methods of deriving the SFR are known to 
differ, in the sense that FIR gives systematically higher 
SFRs compared to Ha. FIR derived SFRs overestimate 
the true current SFR because it is contaminated with 
emission from dust heated by lower mass stars. Ha emis- 
sion is subject to interstellar dust extinction, therefore un- 
derestimating the current SFR. These may lead the reader 
to believe that the discontinuity is not a true one due to 
physical processes. In figure 2, however, the SFR for nor- 
mal galaxies is derived from internal extinction corrected 
Ha luminosity, found by Kewley et al. (2002) to agree with 



FIR based SFR within 10%. The order-of-magnitude off- 
set between the two types of galaxies therefore can not be 
attributed to such observational shorthands. The unac- 
countable offset may be due to differences in the physical 
process of star formation. Elmegreen (1989) has suggested 
that a factor 4 decrease in tg will occur at a shocked sur- 
face of a molecular cloud, and this can account for much of 
the rise in SFR for starbursts, which lies in the "collision" 
regime in figure 2. 

5. Discussion 

5.1. Other parameters 

The simplicity of equation (5) poses a question 
of whether this analysis is not an oversimplification. 
Therefore, we go on to consider other parameters which 
may modify equation (5). The basic assumptions made 
in section 1 is that we have assumed that star formation 
occurs on the free fall timescale tg, but the timescale of 
star formation is still an unresolved issue, and studies of- 
ten give longer timescales than tg. Second, we have as- 
sumed a random motion of the molecular clouds based on 
epicyclic motion of the clouds, while in reality, differential 
rotation in galactic disks is also important to the overall 
motion of clouds. 

The star formation timescale t s f, is equal to tg ~ 10 6 [yr] 
when star formation is dominated by gravitational pertu- 
bations, and this type of formulation t s t ~ tg has been 
adopted in many studies. Other studies argue that t s f 
is longer than tg because magnetic fields provide sup- 
port against gravitational contraction, and star forma- 
tion is impossible without ambipolar diffusion (Tan 2000) 
in which case the star formation timescale is longer by 
an order of magnitude (McKce & Holliman 1999). Still, 
some observational studies suggest that star formation oc- 
curs on scales of ~ 10 6 [yr], based on the observed offset 
between CO and Ha arm in spirals (Egusa et al. 2004). 
This issue is not resolved, and we adopt t s f ~ tg for sim- 
plicity. However, if t s f > tg by an order or more, colli- 
sional star formation will dominate at a much lower ob- 
served density of about E cr ; t = 1.8 — log a. The differ- 
ential rotation of galactic disks causes molecular clouds 
to collide as in the case of random motion. Gammie et 
al. (1991) find a = 5.1 [km • s^ 1 ] for cloud encounters due 
differential rotation, whereas random velocity dispersion 
of GMCs are ~ 7 [km -a -1 ] (See Stark & Brand 1989), 
or ~ 4 [km-s -1 ] (Liszt & Burton 1981; Clemens 1985). 
Therefore, we can say that the value of cloud collision 
velocity, a, is comparable even with the assumption of 
differential rotation. 

5.2. Mergers and interacting galaxies 

The argument of threshold surface density proposed in 
this Letter cannot specify what kind of change will ap- 
pear in the Esfr — S relation such as that seen in figure 
2. The discontinuity may well be a result of transition 
from spontaneous to collisional star formation, but our 
results do not predict or necessitate a discontinuity. A 
sudden rise in SFR as collision dominates, however, is 
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physically intriguing and easy to understand. In a colli- 
sion between two GMCs, the shocked region between the 
clouds arc compressed so that the density abruptly be- 
comes higher in this layer compared to other parts of the 
GMC. This will shorten tg, and consequently the SFR is 
higher (Elmegreen (1989) finds that tg is shortened by a 
factor of 0.25): this should appear as a discontinuity in 
the Schmidt law. 

In this respect, it is reasonable to expect that merging 
and interacting galaxies form stars mainly by cloud-cloud 
collision, and therefore exhibit higher SFR than galaxies 
with the same molecular gas density. In these systems, 
the critical density S C rit is far lower than normal spirals, 
because GMCs in colliding galaxies will have far higher 
collision velocity a. In a collision between two spirals, 
the typical rotation velocity of the galaxies will become 
representative of a, with about ~ 200[km • s" 1 ], depending 
on the sense of rotation of the colliding spirals. In this 
case, logS cr it can be to 1, the typical global surface 
molecular mass density of spirals, so that most of the star 
formation can be attributed to cloud-cloud collision. 

6. Conclusion 

We have constructed a simple model to consider the 
threshold between star formation dominated by spon- 
taneous gravitational collapse and cloud-cloud collision. 
Although both cloud-cloud collision and gravitational col- 
lapse processes have for decades been thought to cause 
star formation, the threshold of these two processes has 
not attracted much attention. Most studies of colli- 
sional star formation have focused on qualitative aspects. 
Numerical studies such as that by Mazzei (1987) have 
pointed out a threshold, but with the contradictory result 
that collisions between clouds result in lower star forma- 
tion rates. 

Our model, by assuming that star formation is a bi- 
modal process dominated either by spontaneous gravi- 
tational collapse or cloud-cloud collision, docs not take 
into account the possiblity of star formation induced by 
expanding SNR shells (self-propagating star formation), 
or complicated structures in temperatures or chemistry. 
Simple as our model is, it succeeds in explaining the ap- 
parent discontinuity in the Schmidt law. Moreoever, our 
formulation coupled with the empirical implication that 
most clouds are virialized and thus GMC size and mean 
density are anticorrelatcd, leads to the conclusion that the 
threshold surface mass density log£ cr j t is practically inde- 
pendent of cloud parameters (size, density, mass), but that 
it depends on the mean velocity dispersion a of clouds. 

Further tests of our formulation can be expected from 
observations of mergers. In merging spiral pairs whose 
rotation are oriented in the reverse direction, the relative 
velocity of the clouds are increased so that according to 
equation (7), collisional star formation starts at a very 
low surface mass density. However, if the rotation is 
oriented in the same direction, a will not increase as much 
as in the reverse orientation. Thus, the threshold log£ C rit 
starts at a higher density. Assuming that collisional star 



formation results in a discontinuity in the SFR relation 
with surface gas density, we can compare mergers with 
these two orientations with average normal galaxies 
which we can assume to be dominated by a spontaneous 
process. Then, we should see that mergers with reverse 
rotation orientation will start its discontinuity in the 
Schmidt law at a lower surface mass density than mergers 
with ordered rotation. 
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